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SUMMARY 


The incompressible rotational flow equations are used to obtain 
solutions for the flow behind axisyrametric shock waves with conic longi- 
tudinal sections . The nonlinear part of the term due to rotation is 
retained in the analysis. Numerical results for standoff distance and 
stagnation point velocity gradient are presented for the case in which 
the shock wave is a paraboloid, a sphere, or an oblate or prolate ellip- 
soid. A similarity parameter is proposed which correlates approximately 
the flow behind geometrically similar shock waves at different free-stream 
conditions . 


INTRODUCTION 


Recently, high-speed computers have been used to great advantage in 
the solution of the supersonic blunt-body problem. The existing methods 
employ either the direct approach (i.e., given a body, the shock wave is 
calculated), reference 1, or the indirect approach (i.e., given a shock 
wave, the body is computed), references 2, 3> and 5- The predictions 
of these theories have been compared with experiment, and results for 
standoff distaftice and body shape (or shock-wave shape) are in good agree- 
ment with the observed values. These methods do, however, suffer from the 
disadvantage that a new calculation must be made for every free-stream 
condition. It therefore seems desirable to obtain approximate analytical 
solutions to the problem which would provide simpler (though less accurate) 
results . These results might be useful in their own right and might also 
aid in the interpretation of the more exact numerical solutions . 

The density ratio across a shock wave becomes very large at high 
supersonic Mach numbers and for ratios of specific heats near 1 (cp/c Y = 

7 si). The subsequent deceleration of the flow to zero velocity at the 
stagnation point can, therefore, change the density only slightly. This 
fact has led several investigators to solve the inviscid incompressible 
rotational flow equations, imposing the exact boundary conditions behind 
the shock wave. Lighthill (ref. 6) has used this technique in the strong 
shock approximation to obtain a solution for the flow behind a spherical 
shock wave. Hida (ref. 7 ) earlier obtained a solution for the flow near 
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the axis. He did, however, consider all supersonic Mach numbers. 

Whitham (ref. 8) , Hayes (ref. 9 ) , and Hida have obtained the correspond- 
ing solution for flow "behind a circular shock wave in two dimensions. 

When the shock wave is a surface of variable curvature, the vorticity 
introduced by the shock wave destroys the linearity of the equations of 
motion. Vinokur, in his application of the method to the study of bodies 
of variable curvature (refs. 10 and 11 ) , has neglected this nonlinearity 
in the term due to rotation. In the present work this nonlinear term is 
retained and a solution is found for the flow behind any axisymmetric 
shock with a conic longitudinal section. For bodies of bluntness zero or 
greater (paraboloid, sphere, or an oblate or prolate ellipsoid) numerical 
results for shock-wave standoff distance, body shape, and stagnation point 
velocity gradient are obtained and comparisons are made with the numerical 
calculations of Van Dyke and Gordon (ref. 4) and Mangier (ref. 5) • Finally 
a parameter is proposed which reduces approximately any blunt body flow at 
a free-stream Mach number greater than 3 "to an equivalent flow at infinite 
free-stream Mach number with a different ratio of specific heats . 

The author is indebted to Dr. Max. A. Heaslet of the Ames Research 
Center for suggesting this problem and for his valuable advice, criticism, 
and encouragement during the course of the research. 
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from body axis 
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distance along body surface 

velocity 
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curvilinear coordinate orthogonal to | 
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density 

azimuthal angle 
stream function 
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value associated with body 
value associated with shock 
value associated with stagnation point 
component in the i) direction 
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component in the | direction 
component In the cp direction 
vector quantity 


FUNDAMENTAL EQUATIONS IN CURVILINEAR COORDINATES 


The constant density approximation as used by Lighthill leads to an 
indirect method in that one begins with a known axisymmetric analytic 
shock wave and investigates the flow behind it. The body shape which pro- 
duced the wave is then calculated. That procedure is followed here. The 
coordinate system is chosen so that the shock wave is a coordinate surface, 
T) = lj cp Is the azimuthal angle; and | is the third member of the 
orthogonal coordinate system. The continuity equation for incompressible 
axisymmetric flow is then: 


div V = 


1 




-3 I 




= 0 


(1) 


where hg, h^, and h^ are the local scale factors of the coordinate 
system. 

A stream function t( | > tj) is defined in such a manner that it satis- 
fies this equation identically. 

Si . . 

S\|r 

It is assumed that vorticity is generated only at the shock wave. It 
then follows that for incompressible axisymmetric flow 

| - ?(t) (3) 

where oj is the vorticity, r is the distance from the axis of symmetry, 
and ?(i) is some function of i alone. This means that if w/r is 
known at some point on a streamtube, it is known everywhere on that stream- 
tube. Thus if f(f) can be evaluated somewhere in the flow field, equa- 
tion (3) can be used to relate (upon substitution for to) the velocity 
components and hence the stream function to a known right side. Lighthill 
has given a result for the vorticity jump across any shock wave (at 

=t») as a function of the shock curvature, density ratio, and component 
of velocity tangential to the shock wave. Hayes (ref. 12) has shown that 
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this result is independent of free-stream Mach n umb er . Therefore f(\|r) 
can he regarded as a known function once the geometry of the shock is 
specified. When the shock wave is axisymmetric , f(\|r) has only a cp 
component. When the indicated substitutions are made in equation ( 3 ), 
the following equation results 

e Cp / h | 

rh^i . d| Vicph^ dr) ytyh^ dr/- 



The Equations in Van Dyke’s Confocal Coordinate System 


Van Dyke (ref. 3) introduced a unified axisymmetric orthogonal 
coordinate system that passes continuously from hyperboloidal to oblate 
ellipsoidal . Thus every conic coordinate system can be obtained from it 
as a special case. It is obtained by setting 


x 

Rs 



^(1 - B s | 2 ) (l - B s + B gl f) 


(5a) 


= I I in (5b) 

Rs 

where Bg is the bluntness of the shock wave. The bluntness , B, of a 
conic can be related to the eccentricity of that conic in the following 
way: 

For hyperbolas and prolate ellipses: 


B = 1 - e 2 


For oblate ellipses: 



(6a) 

(6b) 


For both oblate and prolate ellipses , B is equal to the square of the 
ratio of the r axis to the x axis. A shock wave of bluntness B a 
described by 
oloidal when 


is 


T) = 1. The shock wave is hyperboloidal when B s < 0, parab- 
B s = 0 , prolate ellipsoidal when 0 < B s < 1, spherical when 


•'•Actually for ellipsoids, B s > 0, equations ( 5 ) give only the left 
half of the figure. This is the only part which may form a shock. 
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Sketch (a) 



Sketch (h) 

B s ss 1, and oblate ellipsoidal when 1 < B s * This is illustrated in 
sketch (a) • The azimuthal angle is cp and the third member of the coordi- 
nate system is The coordinate system is shown in sketch (b) . In the 

subsequent analysis, the shock-wave radius is taken to be unity. The line 
element Is given by 


ds 2 = 


L 1 - B s | 2 


(-1 - B s )f 2 + rf 

(1 - B s ) + B s t] 2 


dr f + ^T^dcfi 2 


( 7 ) 


Thus 
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(1 - B s )g 2 + rf 

1 - B,! 2 


h n = 


l ( 1 - B s )| 2 + f 

1 - B s + Bgtf 


(8a) 

(8b) 


h(p = in (8c) 

The function f(\|r) will now "be evaluated. For an axisymmetric shock 
wave in a uniform flow, Lighthill’s result for the vorticity after the 
shock is 


GO 


'<P 


(K - l)‘ 
K 


V t k 


(9) 


where V-^ is the component of free-stream velocity tangential to the 
shock wave in the x, r plane; k is the curvature of the shock wave in 
the x, r plane; and K is p/ P m ‘ Using the definition of curvature, 
the fact that the shock wave is given by equations (5) with r) set equal 
to unity, and equation (9) we obtain 


f(lr) 


-* (K - l ) 2 
~ e cp k 


V„ 


[1 + (1 - Bjr 2 ? 


At the shock wave 


( 10 ) 





Therefore at the shock wave 




v„ 


1 + (1 - B B ) p * 


( 11 ) 


However, since f(\|r) must remain fixed on streamlines, this expression is 
applicable to the entire flow field after the shock. Equation (4) then 
becomes, on substitution of equations (8) and (ll): 
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f 1 - B s | 2 > 

" /2 3x" 

+ A 

A. | 

Si - Bs + Bsn^ 

ax’] 


-IB 

Va - Bg + BgTj 2 / 

3|J 

8r] 

M 

V 1 - Bgl 2 J 

3lJ 


Mt][( 1 - B s )g 2 + rf j 

[(1 - B s | 2 )(l - B s +B s r f)] l/2 [1 + (l - B S )X] 2 



where 


X = 2(K - l) 2 

The boundary conditions, which express the continuity of mass and the 
tangential component of velocity across the shock, are: 

X (I, 1) = I 2 (13a) 

^ (£, 1) “ 2K| 2 (13b) 

OT] 

A solution of equation (12) satisfying the boundary conditions (eqs. (13)) 
would give the details of the flow behind any shock wave of conic longi- 
tudinal section (subject, of course, to the inexactness of the constant 
density approximation) . This equation is, however, nonlinear because of 
its right side and an exact solution to it would be very difficult to 
obtain. It is reasonable to seek a series solution near the stagnation 
point because here the constant density approximation is most faithful. 


Series Solution near the Stagnation Point 


To obtain a series solution we set 


*-> l 2 n fnU) 


(14) 


where the odd powers of | in the expansion can be suppressed by virtue 
of the geometric symmetry about the x axis. This expression is then 
substituted into equation (12) and coefficients of like powers of £ are 
set equal to zero. The following infinite set of ordinary differential 
equations is then obtained: 

(1 - B S + Brf > ± (i li) + b e i£i + i (as. - avi) - A, 3 (l 5 a) 

(1 - B S + B b t,=) A (i fe) + Bs 1 (2^3 - 12B s f E ) 

= X(1 - B s )(r 1 - 2^) (15b) 
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(1 - B s + B s rf ) 


_a_ 

dr) 


if, 


df. 


+ Bs ~S) + t (Wf4 " 30B s f 3 ) 


- *(1 - B s )[3(l - B s )rff!2 - 2(1 - B s ) T]fi - 2rj3f 2 ] (15c) 

etc. On substitution of (l4) into equations (13) the boundary conditions 
become: 

fi(l) - 1 (16a) 


^(l) = 2K (l6b) 

dq 

f ± (i) = ~ (l) « o , i > l (16c) 


A coupling between this set of ordinary differential equations is 
seen to exist, except for B s = 1. In this latter case the right sides 
of all the above equations except the first become zero and this in 
conjunction with the boundary conditions implies that 2 


f i 3 0 , i > 1 

When B g f 1 it would seem, however, that a knowledge of f 2 is necessary 
to calculate f 1; etc. A-power series expansion about q = 1 is, however, 
possible. Since f 2 (l) and df 2 /dq(l) are both known from the boundary 
conditions (are in fact zero), d 2 fi/dq 2 (l) and d 3 f;]_/dq 3 (l) can both be 
evaluated from the first equation, and d 2 f 2 /dq 2 (l) and d 3 f 2 /dq 3 (l) can 
then be evaluated from the second equation. This process can be carried 
on indefinitely. The following series expansion for X is then obtained. 

x = i s fi(n) + £ 4 f 2 (ii) + £ 6 f 3 (n) + ^(q) + i 10 f 5 (q) + • . . (17) 

where 


2 Vinokur*s assumption that the nonlinearity in the term due to rota- 

tion is unimportant also leads to this result and he therefore needs to 

integrate only equation (l5a) . 
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fiU) - 




f3(n) 


1 - 2K(l - rj) + [-(K - 1)B S + (K 2 - K + l) ] (l - T]) 2 + -(K - l)B £ 

(1 - n) 3 +i (K - 1) 


+ | (K - 1)(3K - 2)B S - | (K - 1)' 


~ 5BS" 


+ (5K - 3)B S 2 - | (3^K - 35)B S + ^ (K - 1) 


(1 - nr 


+ I (K -1) 


1 B s 4 + | (7K + 3 )b s 3 - i (280K - 283 )B S 2 


+ JL (32K 2 + 2 k 5 K - 275)B b - i (k - l) (i6k + 47) 


(1 - n) 5 


+ (K - 1) 


_ 21 B s 5 + I (9K - 2)B S 4 - \ (189K - 186 )B S 3 + ± (l6 OK 2 

2 b 6 30 


+ 987K - 1144) B s 2 + i (- 32K 3 - 736K 2 - 1236 k + 2011 )B S 


+ — (K - 1) (16K 2 + i 44 k + 251) 

45 


(1 - n) 6 + • • 


(18a) 


(k-i) 2 ^(Bs- i)(i-n ) 2 + 


B s 2 - | (2K + 5)Bs + ~ (4 k+ 7) 


(l - 7)' 


+ — [15B S 3 - (24K + 6l)B s 2 + (4K 2 + 56K + 101)B S 

- (4 k 2 + 32K + 55) ] (l - n) 4 + — [105 B S 4 - 5(36k + 95)Bc 

60 

+ (4 ok 2 + 612K + 983)b s 2 - (112K 2 + 700K + 985 )b s 


+ 4 (i 8 k 2 + 67K + 93) ] (1 - n) 5 r + • • • 


(18b) 


(K - 1) 2 4 (B s - 1) 2 (1 - n ) 2 + I [3B S 2 - (8k + 15)B S 2 + (16K + 21)B S 


- (8k + 9)](1 - n) 3 + 3^ [ 15 B S 4 - 8 (6k + ll) B s 3 + 2(i8k 2 + 84 k 
+ 115)B 2 - 4 (i 8 k 2 + 48 k + 64 )b s + 9(4 k 2 + 8k + ll)](l - n)^ 


(18c) 
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f 4 (ri) = (K - 1) 



- 1) 3 (1 - T ]) 2 + | [ 3 B S 4 - 


4(3K + 5)B s s + 6(6 k + 7)B S 2 


- 36(K + 1)B S + (12K + n)](l - 



(l8d) 


fs(B) = (K - 1) 2 (B s - 1) 4 (1 - rf) + . . . (l8e) 

The question of convergence of the double series now arises. Since 
no general expression was found for the terms in the series, it was not 
possible to investigate the convergence directly. It can he verified by 
direct substitution, however, that for | B s | <1, the terms decrease 
sufficiently rapidly that the series appears to be of practical value 
for computation. For | B s | >1 the terms decrease so slowly that the 
series is of little value. A transformed version of this series will 
therefore be given which does converge sufficiently rapidly for all 
positive B s . 


Transformation of the Series 


The transformations 


I = V 1 - B s 


(19a) 


T) = 4 BgT) 2 + 1 - B e 


f = — f 
x n jam 


(19b) 

(19c) 


when applied to equations (15) give the following result: 

A ,_ s 

_ •— + 1 o±‘o — yr-. ) — = ■ 

dip 


sr + <&* - itr (V - 1 * B S ) 


(20a) 


+(24f 3 -12f 2 ) = 


q 2 - 1 + B s Be 


[1-B S - 2(1 -B s )0f - 1+B s )fi] 


( 20 b) 


+ (48f 4 - 30f 3 ) - X 

dq 2 


T - 1 + Bg 


_A 

B s 


■4 [-2(l-B s )fi -2(l-B s )(Tl 2 -l + B s )f 2 

2 / —r*. „ „ \ TT 2 • 


+ 3(1 ~ B s ) (rp-l + Bjfi ] 


(20c) 


The boundary conditions become, using (19) and (16) 
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fi(l) 


_1_ 

B s 


df i fi) 2K 
drj 1 “ B S 2 

fi(l) a ^ (l) « 0 
dn 


i > i 


(21a) 

(21b) 

(21c) 


The following double series for X can now be obtained either by 
transformation of equations ( 17 ) and (l8) or by straightforward solution 
of equations (20) in the manner previously indicated. 


x= (1 - i 2 )f 1 (r 1 ) + (1 - ? 2 ) 2 f 2 (?,) + (1 - l 2 ) s f 3 0 i)(i - I 2 ) 4 f 4 (n) 

+ (1 - I 2 ) 5 f 5 (n) + . . . (22) 

where 


fi(^) “ B s _1 - 2K(1 - ^)B S - 2 + B s " 2 [1 + (K - l) S B s -1 ] (l - ^) 2 

- | (K - 1)B S “ 3 [1 + (K - l)B s -1 ] (l - ^) 3 

- | (K - 1)B S - 4 [(K + 1) - 2(K - l)B s - 1 ](l - Tj) 4 
+ — (k - i)b s -4 [i + 4k 2 - 5K - 5)b b " 1 

- 2(K - 1)(2K - l)B„ -1 ](l - Ti) 5 + — (K - l)B q -5 [3(2K + l) 

4-5 

- (2K 3 - 14k 2 + 37K - l)B s -1 + 2(K - l)(K 2 - 6 k + 11)B s “ 2 ](1 - ^.) 6 


+ . . . (23a) 

f 2 (q) = (K - 1) 2 |(B s - 1)B s ~ 4 (1 - ^) 2 - I B s " 4 [(1 - Bg" 1 ) (K + l)](l - ^) 3 

- | Bg^El - (K 2 + 2K + 8) + (K 2 + 2K + 7)B s “ 2 ](l - ^) 4 

+ ~ B g ~ 5 [ (i4k + 1) - (4k 2 + 30K + 5)B S _1 
+ (4K 2 + 6K + 4)B S " 2 ](1 - ri) 5 j + . . . 


(23b) 
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fa(Ti) 


= (K - 1) 2 |(1 - B s 1 )B S “ 3 (1 - T ]) 2 - | (1 - b s - 1 )b s -*(4k + 3) (1 - n) 3 


- i B„ 4 [1 - (9K 2 + 6K + 15) - (9K 2 + 6 k + 17)b s “ 1 


+ (18K 2 + 12K + 31)B S ~ 2 ](1 - T)) 4 ^ + . . . 


(23c) 


f 4 (n) « (k - i) ; 


(1 - B s X )Bs 3 (1 - ij) 2 


fs(tl) - (K -1)' 


~ (1 " B S - 1 )B S - 4 (3K + 2)(1 - ^) 3 J + . . . 

n) 2 ] + • • • 


(1 - b s -1 ) 4 b s - 3 (i 


(23d) 

(23e) 


The following computation was undertaken to obtain information about 
the convergence of the series in equation ( 22 ) , The first of equations 
(20) was integrated numerically assuming f 2 =_0. The first two equations 
were then integrated simultaneously assuming f 3 = 0. This process was 
continued through the simultaneous integration of the first five equations. 
The results were then compared with each other and with numerical values 
given by the series. A typical result is shown in table I. 

It can be seen from table I that f^ is only slightly affected by 
the addition of the other f's. It is, therefore, reasonable to suppose 
that the results from the simultaneous integration of the first five of 
equations (20) are for all practical purposes the exact solution for f^. 
The series solution is then seen to be extremely close to the exact 
solution. 


CALCULATION OF FLOW QUANTITIES 
Standoff Distance and Body Bluntness 


The stream function has been so chosen that the axis of symmetry is 
the zero streamline. When this streamline reaches the body it wets the 
entire surface. The body can then be located by examining the zeros of the 
stream function. Near the axis, | = 1, 

x ~ f x (n) 

The body shape near the axis is therefore given by a zero of fi, and 
since f x is a function of r| alone, the body is necessarily a confocal 
ellipse. 
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Away from the stagnation point, 1 - g is not_very much less than 
1, and the body surface is not given by a zero of f^. The body coordi- 
nates are solutions of the equation: 

fiCq) + (l - 5^ ) ( ^)) +•••** 0 

The magnitudes of f 2 , f 3 , . . . have been investigated both from the 
series obtained and by numerical integration of the differential equations. 
It was found that these functions are sufficiently small that, in the 
region where the initial assumptions can be expected to hold, they are 
negligible. The body is considered, therefore, to be a conic. Following 
references 3 and 5 bodies obtained here are characterized by a bluntness 
Bg. The bluntness of an ellipse defined by rj “ rL, is (from eqs. (5), 

(6), and (19b)) 




T| b s — 1 Hh Bg 


(24) 


This quantity together with the body nose radius completely specified the 
body shape. The body radius Rd when referred to the shock-wave radius 
can be expressed as (from eqs. "(5) and (19b) and the definition of radius 
of curvature) 

% - 1 (25) 
Rs 


For very large K, an approximate zero of f x can be found (from eq. 
(23a) ) . When this value of 11 is substituted in equation (ja) and in 
equation (24) , the following relations result: 


J>_ = 1 _ /§ 1 £ ~ % 

RB " K si 3 K 372 K 2 


( 26 ) 


R B_ 1 B B“ 2 [ 8 ^ ~ 2 

R s K J 3 K 3/2 


(27) 


Equation (26) is in agreement with the result given by Hayes 3 (ref. 9) for 

*D — 1 

JDt5 X * 


If the shock-wave shape is held constant and K varied, the body 
shape will, in general, change. Figures which show standoff distance and 
shock radius for a given body as a function of K are, therefore, cross 


3 Since Hayes presents &/r s , equation (27) must be used to show that 
(26) is identical to Hayes’ result for Bg = 1. 
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Sketch (c) 

plots. Sketch (c) compares standoff* distance calculated "by the present 
theory to that given in reference 5 for a shock of hluntness one half. 

At K = 6(7 = 1.4), the present calculation differs from reference 5 "by 
5-8 percent, whereas at K - 21(7 = l.l), the difference is only 
0.3 percent. Thus it seems that for small 7 the constant density 
approximation is very accurate. Figure 1 gives the ratio of standoff 
distance to "body nose radius as a function of K for a large number of 
bodies. Vinokur’s results are shown for comparison. 

Figure 2 compares body radius obtained from the constant density 
approximation to results of reference 4. It can be seen that the constant 
density approximation predicts a larger body radius than does the exact 
computation . 

Figure 3 gives body bluntness as a function of shock-wave bluntness 
for several 7. Further blunt-body solutions of the type presented in 
reference 4 have been obtained at the Ames Research Center by 
Mr. Franklyn B. Fuller and two of these results are plotted for comparison. 
It can be seen that the exact computation predicts a larger body bluntness 
than does the constant density approximation. This disagreement decreases 
for smaller 7, and for 7 = 1 both methods give the same result. 
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Stagnation Point Velocity Gradient 


The gradient of the velocity at the stagnation point is 

Rs /dv\ = B s 3 f x »(^) 

Voo \ ds / st 2K( T ) 2 + B s - 1) 


(28) 


In the table below, stagnation point velocity gradient is compared with 
that of reference 5 • As can be seen, the agreement is better for small 
7. It is also evident that the constant density approximation gives a 
reasonably accurate estimate of the stagnation point velocity gradient 
even for 7 = 1 A. 


Moo 

7 

Bs 

Rs / 

Voo ' 

'd\A 

^ds/ s -(- 

Percent 

error 

Hef. 8 

Present 

00 

1.4 

i/2 

0-7645 

0-7564 

1.06 

00 

1.2 

i/2 

0.5332 

O.5296 

0.68 

00 

1.1 

1/2 

0.3723 

0.3712 

0.30 


For large K the following result can be derived for stagnation point 
velocity gradient 


^ - /! _±_ + 2 foz + . . . (29) 

T »W st J 3 k 1 '2 9 7575 k 3 ' 2 

Hayes gives the first term of equation (29) in reference 9- The 
functional form of equation (29) is in agreement with the expansion for 
large K of the classical result given by Li and Geiger (ref. 13) with 
the constants differing slightly. Neither of these latter results includes 
the effect of body bluntness . 

Figure 4- is a plot of stagnation-point velocity gradient as a function 
of body bluntness and K. Vinokur’s results are shown for comparison. 


Pressure Distribution on the Body 


The vector momentum equation for incompressible flow can be written 



V X Curl V = grad 


(30) 
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Scalar multiplication of equation (30) by V yields 

^' V (p + f") = ° (31) 

— ► 

where V*V is a derivative in the direction V and is therefore the 
derivative along a streamline. The above equation expresses the fact 
that 

| + 7T- = (32) 

Equation ( 32 ) is the Bernoulli equation for incompressible rotational 
flow. The pressure variation on any streamline in the flow can be calcu- 
lated from it (presuming that H is known at some point on the stream- 
line) . It is evident from symmetry that the streamline which passes 
through the normal shock wets the body surface . Using the strong shock 
relations to evaluate the pressure, density, and velocity behind the 
normal shock one obtains 


H(0) = 

K 
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The pressure on the body surface is then 


(33) 


P 

(1/2) PcJ 


2 

00 


1 + V 

K 7+1 





and (v/Voo) 2 is then computed from the stream function. Taking account 
of the facts that on the body surface fi(t) is nearly zero, and that near 
the axis (l - I 2 ) « 1 this equation becomes 


+ 


4y 


B e 


1 - r 


(l/2)p 00 V 00 2 K 7+1 ^ If + (B s - !)p 


[fi^n )] 2 


(35) 


Figure 5 is a plot of the pressure distributions on a sphere given by the 
constant density approximation and by Van Dyke and Gordon (ref. 4-) . 


APPLICATION TO FLOW AT FINITE FREE-STREAM MACH NUMBER 


In the preceding analysis it was convenient to assume that the free- 
stream Mach number was infinite because of the fact that as a consequence 
of the shock-wave relations, the density is constant just after the shock 
wave only if the Mach number is infinite. It will be shown, however, that 
near the axis of symmetry, the variation in density is not great at 
moderately high free-stream Mach numbers. 



l8 


The density ratio across a shock wave satisfies the equation 


= LiLi (36) 

Poo 7 - 1 + (2/yi^sln z a) 


where a is the angle of inclination of the shock wave to the free stream. 
It can he seen from this equation that if attention is restricted to a 
small region containing the normal portion of the shock wave (a £ jt/2) , 
the variation in p/p can "be made as small as desired. Then sin a is 
approximately 1 and we have the following relation for K: 


7+1 

7 - 1 + (2/m/) 


( 37 ) 


It can he concluded from the preceding argument and from the fact 
that Hayes has shown that Ligh thill’s result for the vortlcity after a 
shock holds at finite free-stream Mach number, that the results arrived 
at previously can reasonably he applied to a small region of the flow 
behind a shock wave in a flow with finite free-stream Mach number. The 

flow downstream of the shock wave is seen then to he a function of B s 

and K or, for a given shock wave, K alone. This means that (within the 
limits of applicability of this theory) a flow at any 7 and free-stream 
Mach number has an analog at Moo *= 00 and a different 7. Fuller’s numer- 
ical solutions have been used to test the validity of this parameter. 

Figure 6 shows the result of two calculations, one at Mcc ** 3 > 7 = 1-4 
and the other at Moo = 00 and 7 = 1.7* In both cases 

B s = 0.25 


~ ± a 3.857 

7 - 1 + (2/M*, 2 ) 

It can be seen that the two bodies are very similar in shape, and 
that the shock-wave standoff distances are nearly the same . The pressure 
distributions are almost identical. Figure 7 shows standoff distance for 
finite free-stream Mach number and shock bluntness of one-half plotted 
against (7 + l)/(7 - l) in one case and {(7 + l)/[7 - 1 + (2 /Moq 2 )]} in the 
other. The exact solutions are from references 4 and 5 and Fuller. 

In both cases (i.e., body shape and standoff distance) the correlation 
is better at higher free-stream Mach number. 


CONCLUDING REMARKS 


A solution has been presented for the flow after axisymmetrie shock 
waves with conic longitudinal sections. The flow behind the shock wave 
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has been assumed incompressible but vortical, and the nonlinearity in the 
term due to rotation has been retained. Shock-wave standoff distance, 
body bluntness, and stagnation point velocity gradient have been presented 
for bodies of bluntness greater than zero. Pressure distribution on a 
sphere has been presented also. Stagnation point velocity gradient is in 
very good agreement with that presented in reference 5* The agreement of 
standoff distance with results of references 4 and 5 is reasonably good 
( 5*8 percent) at y - 1.4 and improves with decreasing 7 ( 0.3 percent 
at 7 = l.l). Thus it can be concluded that at high supersonic Mach 
numbers and for 7 near one, compressibility is unimportant, but vorti- 
city is important . 

The quantities mentioned above, standoff distance, stagnation point 
velocity gradient, and body bluntness, agree with those presented by 
Vinokur. It therefore appears that he was justified in neglecting the 
nonlinear part of the term due to rotation. 

A parameter has been proposed which approximately correlates the flow 
after geometrically similar shock waves at different free-stream Mach 
numbers. It is hoped that this parameter may prove useful in cataloging 
numerical blunt-body solutions. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 27 , 1961. 
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TABLE I.- COMPARISON OF RESULTS OF NUMERICAL INTEGRATION 


WITH THOSE OBTAINED FROM POWER SERIES 
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Figure 4.- Stagnation point velocity gradient at ^=00 for several values of bluntne 

parameter, Bg* 
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